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Abstract 

We use the generalized zeta function regularization method to compute the one-loop 
quantum correction to the masses of the TKl and TK2 kinks in a deformation of the 0{N) 
linear sigma model on the real line. 
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1 Introduction 



In this paper we shall apply the generalized zeta function regularization method to unveil the semi- 
classical behaviour of the quantum topological kinks arising as BPS states in multi-component 
scalar (l-l-l)-dimensional field theory. The interest of this investigation lies in the fact that these 
systems live at the heart of the low energy regime of string/M theory. Effective theories such as 
N=l SUSY QCD and/or the Wess-Zumino model encompass A^-component scalar fields and 
have a discrete set of vacuum states. From a one- dimensional perspective one can foresee the 
existence of fc-component topological kinks, k < N, which are the seeds of the BPS domain walls. 
These extended states play such a prominent role in the three-dimensional world that analysis of 
the quantum behaviour of fc-component topological kinks becomes an important issue. 

Here, we shall focus on a one-parametric family of deformations of the 0{N) linear sigma 
model, in (l+l)-dimensional space-time. The model to be addressed forms the bosonic sector of 
a super-symmetric theory with real super-fields and is of the general type of the Wess-Zumino 
model, like those studied in and 0. An important feature common to all these models is the 
non-existence of continuous symmetries in internal space and they can therefore be included in the 
class of systems considered in Reference 0. The authors of these works use continuous phase-shift 
methods to calculate the one-loop correction to the kink energy in full super-symmetric theories. 

Our approach differs in two ways: 1) we restrict ourselves to the bosonic sector and leave the 
fermionic fluctuations for future research; 2) the generalized zeta function regularization procedure 
is applied to deal with the infinite quantities arising in the quantization prescription. This method 
has been used previously in the description of quantum corrections to kink masses for theories 
with a single real scalar field, see but here we shall follow the more elaborated procedure 
developed in p and for one-component systems. Nevertheless, our method is particularly 
suited to models with several real fields because in general, in theories of this type, the spectrum 
of the second order fluctuation operator around kink solutions is only partially (asymptotically) 
known. Since this operator is a Schrodinger operator acting on functions belonging to the Hilbert 
space L^(R) ® C^, one can write an associated heat equation. From the asymptotic expansion 
to the heat function one obtains enough information about the generalized zeta function of the 
second fluctuation operator |]T0| to provide a good approximation to the one-loop correction to 
the kink mass. Asymptotic methods to compute the mass of quantum solitons were first used in 
Reference . 



Within a given range of the deformation parameter, the modified 0(2) linear sigma model 
that we shall consider is the celebrated Montonen-Sarker- TruUinger- Bishop ( MSTB ) model 
0. Over the years, several kinds of kink solitary waves with very noticeable properties have been 
discovered in this system , [Q. The kink moduli space of the A^ > 3 generalizations of the MSTB 
model has been described in ||^ , whereas more recently the stability of the different kinds of kink 



was established in |T^. Therefore, computation of the semi-classical mass of the stable topological 
kinks with one non-null component, TKl , or two non-null components, TK2, is compelling. To 
achieve this task, which is the main goal of this paper, the zeta function regularization method is 
specially appropriate , because of the impossibility of solving the spectral problem of the second 
order fluctuation operator. 

The organization of the paper is as follows: In Section §.2 the general semi-classical formula 
for the mass of quantum solitons, the zeta function regularization procedure, and the zero point 
energy and mass renormalization prescriptions are explained. In Section §.3 we compute the one- 
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loop quantum corrections to the masses of topological kinks with one and two non-null components 
in the MSTB model. In Section §.4 similar formulas are given in the deformation of the linear 
0{N) sigma model, with > 3, that generalizes the MSTB model. Section §.5 offers some 
comments on possible extensions of our results. Finally, in Appendix A the first coefficients of 
the matrix heat kernel expansion are written, whereas, in Appendix B, it is shown that only the 
stable two-component topological kinks saturate the topological bound. 

2 Semi-classical mass formula for k-component quantum 
topological kinks 

A^-component scalar fields are maps from the M}'^ Minkowski space-time to the "internal" 
space: 

N 

. Tn>N 



a=l 

Here, y^, fx = 0, 1, are coordinates in R^'-*^; = is a basis in TR-*^'^, and Ca, a = 1, 2, • • • , A^, 
are ortho-normal vectors in M^, Cq ■ Cf, = 5ab- We shall consider (l+l)-dimensional field theories 
whose classical dynamics is governed by the action 
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S = j d'y^-d,i;d^i;-U{i;)y 

We choose the metric tensor in T^(M^'^) as, g = diag(l, —1) and the Einstein convention will be 
used throughout the paper only for the indices in M^'^. The system of units -identical to that 
chosen in 0- is such that only the speed of light is set to c = 1. 

The classical configuration space C is formed by the static configurations ipiy) - we denote the 
spatial coordinate a.s y^ = y - for which the energy functional 



i4t 



is finite: C = {ip{y) / E{ip) < +oo}. In the Schrodinger picture, the quantum evolution is ruled 
by the Schrodinger functional equation 

zh^^^[i;{y),t = y°] = Hvl>[^(y),t = y\ 
If ^pv is a constant minimum of U, the masses of the fundamental quanta are: v'^Sab ~ ' ^ 



The quantum Hamiltonian operator 
acts on wave functionals '^[^{y),t\ that belong to L'^{C). 




3 



A straightforward generalization of the arguments and definitions of the Section §.2 of |0] - 
following the classical papers |Tl|, |jl2[ - shows that the kink ground state energy at the semi- 
classical limit is: 



= E[M + ^ Tr (P/C)^ + o{h^) = E[^^] + ^ Cp^(-i) + 0^ . 
Here )C is the second variation operator, with spectral equation /C^„(x) = co'^^„(x), and 



(1) 



t<;2>0 



2\s 



is the associated generalized zeta function; P is the projector to the strictly positive part of Spec /C. 
In this case, however, K, is a N x N "matrix" differential operator with "matrix elements" 



ab 



dab + 



dy 



6lpaSlpb 



a,b = 1, 2, 



N. 



2.1 Generalized zeta function regular izat ion method 

We shall regularize Cp/c(~|) by defining the analogous quantity Cpic{s) at some point in the 
complex s-plane where Cpic{s) does not have a pole. Cpic{s) is a meromorphic function of s such 
that its residues and poles can be derived through heat kernel methods, see |10]. If Ki^{y, z; (3) is 
the kernel of the N x N "matrix" heat equation associated to /C, 



d 



—lN + IC]Kjc{y,z;l3)=0 



Kic{y,z;0) = iN^iy - z) , 
the Mellin transformation tells us that 

r(s) Jo 



(2) 



where 



J 



-J + 



dyK/c{y,y;P) 



is the heat function hpic[(3], if /C is positive semi-definite and dimKer/C = j. The "regularized" 
kink energy is at the semi-classical limit: 



(3) 



where /x is a unit of length"^ introduced to render the terms in @ homogeneous from a dimensional 
point of view. The infiniteness of the bare quantum energy shows itself here in the pole that the 
zeta function develops for s = — |. 

The renormalization of lE^(s) will be performed in the same three steps as in 0: 

A. The quantum fluctuations around the vacuum are governed by the Schrodinger operator: 



V, 



ab 



dab + 



dy 



where ipv is a constant minimum of U and Uabi'ipv) — jtjrs^itpv ~ 
variational derivatives of U at ipv The kernel of the heat equation 

d 



vl5nh is the matrix of second 



d(5 



provides the heat function hy{(3), 

hviP) 



Tr e 



-/3V 



N 

E 

a=l 



and, through the Mellin transformation, we obtain 

1 



Cv{s) 



Tis) 



df3 f]'-^ Tie- f^^. 



The regularized kink energy measured with respect to the regularized vacuum energy is thus 



E 



Em + ^f^''^' iCpds) - Cv{s)] + o{h' 



B. If we now pass to the physical limit E = lim^^_i E (s), we still obtain an infinite result. 
The reason for this is that the physical parameters of the theory have not been renormalized. 
It is well known that in (l+l)-dimensional scalar field theory normal ordering takes care of all 
renormalizations in the system: the only ultraviolet divergences that occur in perturbation theory 
come from graphs that contain a closed loop consisting of a single internal line, . From Wick's 
theorem, adapted to contractions of two fields at the same point in space-time, we see that normal 
ordering adds to the Hamiltonian the mass renormalization counter-terms 



H{5m^) = 

up to one-loop order. To regularize 

6 ml 



2 y ^s*:^*:^^''*'^ 



dk 



5n 



4tt 



we first put the system in a ID box of length L so that Sm^^^ = ^Cvaai^)^ if the constant 
eigen-function of Vaa is not included in Cv„„. Then, we again use the zeta function regularization 
method and define: (5m^^(s) = —j^^^y-fJ^'^^'^^Cvaai^'^ Note that Sm^^ = lim^^_ i 5m^^(s). 
The criterion behind this regularization prescription is the vanishing tadpole condition, which is 
shown in Appendix B of Reference to be equivalent to the heat kernel subtraction scheme. 
The one-loop correction to the kink energy due to H{5m^{s)) is thus 



hm — /i 



H{5m\s)) ^k) - {^v\H{5m\s)) 

N 



2s+ 



— lim — /i 

L— >oo zIj 



2s+ir(5 + i 



s + l) 



Tis) 



^ 2 a=l 
N 

a=l 



5^U 




6^U 













dyVaaiv) 



(4) 
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because the expectation values of normal ordered operators in coherent states are the correspond- 
ing c-number-valued functions. 

C. Finally, the renormalized kink energy is : 

= E\'^k] + lim [Ai£^(s) + A2£^(s)] + o{h^) . (5) 

2.2 Asymptotic approximation to semi-classical kink masses 

In general it is very difficult to compute /ipk:[/3] exactly. In such a case, we shall make use 
of the asymptotic expansion of /ip/c[/5], which is well defined if < /5 < 1. In order to use 
the asymptotic expansion of the generalized zeta function of the /C operator to compute the 
semi-classical expansion of the corresponding quantum kink mass, it is convenient to use non- 
dimensional variables. We define non-dimensional space-time coordinates = mdU^ and field 
amplitudes (j){x^) = Cdip{y^), where and Cd are constants with dimensions [nid] = and 
[cd] = M-^L-h. Also, we write [7(0) = ^U{'ijj). 

d 

The action and the energy can now be written in terms of their non-dimensional counterparts: 

The important point is that the Hessians at the vacuum and kink configurations can always be 
written as the N x N matrix differential operators V = — and K, = — ^IC, given by 



X 



ICab = '^^'^b + vl6ab - Vab{x) 



where t:^]^^ = and ^^^1^^^ = ^a^ab - Vab{x). Therefore, 

1 . . . . , . 1 
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Cv(s) = — Cv(s) , Ck{s) = — C^(s). 



The asymptotic expansion of N x N matrix heat kernels is given in Reference [|T5|. Nevertheless, 
we shall sketch this procedure in order to adapt it to our computational needs. We thus write the 
heat kernel in the form 

K;^{x, x'; f3) = A{x, x'; (3)K^{x, x'; 13) 



where 

[K^]ab{x,x'](3) = -j==exp 



(x — x'Y'^ 



4/3 



exp{-Pvl)6ab 
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is the solution of the heat kernel equation for the V operator with the initial condition x'; 0) 

6abS{x—x'). The matrix elements of the NxN matrix- valued function A{x, x'; j3) satisfy the system 
of N"^ coupled PDE: 



N N 

EE 

c=l d=l 



d 92 



d[3 dx"^ 



(6) 



whereas A{x,x'; 0) = l^v is the NxN identity matrix. 

For /5 < 1, we solve (^) by means of an asymptotic (high-temperature) expansion: A{x, x'; [3) 
Yl'^=o^ri{x,x')P"'. In this regime the heat function is given by: 



N 



N oo 



Tre"^^ = ^ / dx[Kj 

a=l -^-^^ 

N oo _g^2 

- EE' 

a=l n=0 



a=l n=0 * ' 



oo 
-oo 



[a„]aa(£)r. 



The coefficients [a„]a6(a;, x') can be found by means of an iterative procedure that starts from 
[O'o]ab{x,x') = 6ab and becomes more and more involved with larger and larger N. In Appendix 
A, this procedure is explained and the explicit expressions for the lower orders are shown. 

The use of the power expansion of hpj^lP] = Tre~^^^ in the formula for the quantum kink 
mass is developed in three steps: 

1. First, we write the generalized zeta function of V in the form: 



1 mgL 
T(s) Ja^ 



N 



a=l 



J2 / d(3(3'-h-f''"+B^{s) 



with 



N 



Bv(s) 



a=l " 



2s-l 



r r 1-2 



a=l 



Tis) 



r[s,t>^] and 7[s — being respectively the upper and lower incomplete gamma functions, see 



T6| . It follows that Cvi^) ^ meromorphic function of s with poles at the poles of 7[s — 



which occur when s — | is a negative integer or zero. B^{s), however, is an entire function of s. 

2. Second, from the asymptotic expansion of hpj^[P] we estimate the generalized zeta function 
of ^: 



r(.) 



Jo V 47r ~- y--^ Jo 



a=l n<no 



+ B 



PK\ 



N 



+ 



sT{s) r(s)v^ 



XI 5Z Maa(/C) 



- 7[s + n- i, 



2' a\ 



a=l n<no 



_2(s+n-i) 
Va 



where 



^ N oo r . 



2 ' ^al 



a=l n>no 



Va 
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is holomorphic for Res > —no + ^, whereas 



B 



PK\ 



1 ^ 

J7\ XI 



r(s) 



a=l 



d/?Tr[e- 



is an entire function of s. The values of s where s + n — | is a negative integer or zero are the 
poles of Cp£(s) because the poles of 7(5 + n — \iv'^ lie at these points in the complex s-plane. 
Renormalization of the zero point energy requires the subtraction of Cvi^) from Cpici^)- We 



find: 



Cpk{s) - Cv{s) 



T s) 



a=l n=l ^ Va 



and the error in this approximation with respect to the exact result to Ai£^ is: 

erron = ^[-^ h^,^jc{-\) + Bpd-\) ' Bv{-\)]. 

Note that the subtraction of Cv{s) exactly cancels the contribution of ao(^) and hence the diver- 
gence arising ats = — |,n = 0. 

3. Third, Ai£^ now reads : 



h 



2^ - - ■ 2 

N no-1 



TV 



a=l 



[ai]aa(/C)7[g+ jt^g 



47rr s 



+ 



EE 



l2n-2 



- a=l n=2 v/4^r(-i) e 

The surplus in energy due to the mass renormalization counter-term is: 



AT 



lim > 



a=l 



2L 



2s+l 



r(.) 



■Ck.(s+1) + o(/i2) 



2V47r 



a=l 



2s+l 



and the deviation from the exact result is : 



N 



errors = Yirn^ — Xl[«i]aa(^)^Ka(i 



a=l 



Therefore, 



2v^ 



TV no-1 



^ a=l n=2 « 



+ 



-F _^"!i y^[ai]aa(/C) lim^ 



2V47r 



0=1 
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The contributions proportional to [ai]aa{l(^) of the poles at s 
exactly cancel. 



i in Ai£:-^(s) and Aae^fs) 



We thus obtain the very compact formula: 



J 



AMk = hma [Ao + D„J < 



D 



no 



2v^ 

N rto-l 

-EE 

a=l n=2 



(7) 



6T[ 



a 



In summ: there are only two contributions to semi-classical kink masses obtained by means of 
the asymptotic method: 1) hnidAo is due to the subtraction of the zero modes; 2) hrndDn^ comes 
from the partial sum of the asymptotic series up to the uq order. We stress that the merit of the 
asymptotic method lies in the fact that there is no need to solve the spectral problem of /C: all 
the information is encoded in the potential V{x). 



3 The MSTB model: a deformed 0(2) linear sigma system 

Let us consider now a = 2 model determined by the potential energy density, [^: 

u[^] = \{^i^-'^? + ^i'l- (8) 

The system is a generalization of the A(0'^)2 model to a = 2 scalar field, the 0(2)-linear sigma 
model, although it has been deformed by a quadratic term in ■ip2 in order to avoid the Goldstone 
boson. The deformation parameter a has dimensions of inverse length, [a] = L^^, and the choice 
of non-dimensional variables in the so called MSTB model comes from the choice of = — , 
TTT-d = ^ and = rn?a'^ : 

(T^ is the non-dimensional parameter that measures the deformation and we shall also use the 
related parameter o"^ = 1 — o"^. f/(0) is not invariant under the 0(2)-matrix group if cr^ is not 
zero. The "internal" symmetry group is the Z2 x Z2 group of discrete transformations: 

The vacuum classical configurations are: 

(pv{x,t) = ±ei , iJv{y,y^) = 

Therefore, if we denote by 

the Z2- subgroup that leaves 65 invariant, the vacuum orbit and the vacuum moduli space are 
respectively: M = -%)- = Z2, M = ^ = point. 

There are two kinds of topological kinks which are thus loop kinks and candidates for being 
stable solitary waves of the system: 
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Topological kinks with one non-null component: 



(PTKi{x,t) = (tanhxjei , ipTKi{y,y ) = ^(tanli^)ei 



• Topological kinks with two non-null components: 

<pTK2{x,t) = [(tanhcrxjei ± (T(sechcra;)e2], 

^TK2(2/,l/°) = --^[m(tanh^)ei ± (m^ - a2)(sech^)e2]. 
Note that there are two two-component topological kinks , which only exist if < cr^ < 1. 

The kink and vacuum solutions have classical energies: E[t/vpKi] = f:^) E[^tk2] = 2a{l — 
and Eltpv] = 0. Thus, E['^tki] > E['j/Stk2] and we shall see that the TKl kink, as a 
quantum state, is unstable if cr^ < 1. The lower bound in energy in the topological sector of 
the configuration space is attained, however, by the TK2 kink. This last statement is proved in 
Appendix B. 

The Hessian operator for the vacuum solution is: 

£,+2m^ \ m^/-^ + 4 A m^^ 



and hence the masses of the fundamental quanta are vf = 2m — and V2 — ^ — ^^2, in 

(2) 

such a way that the H} symmetry is spontaneously broken. 
The Hessian operators for the topological kinks read: 

• TKl, 

_ / -4. + 4 6 

)C — { ^ cosh^ X 

2 

TTl — — 

and /C = — ^- Note that /C has a negative eigenvalue if cr^ < 1, as it should be for a 
unstable solution. 

• TK2, 

'17 _ I ax'^ cosh cx cosniTi; 

I 4^tanh£x S I ^2 I 2(2-3a^) 

\ coshfTa; dx^ cosh^ cx 

2 

m - 



and V, — — H. 
2 

The corresponding generalized zeta functions satisfy 
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3.1 The quantum TKl kink: exact computation of the semi-classical 
mass 



• Generalized zeta function of V: 

Acting on the L^{R) (g) Hilbert space, we have that 

SpecV = {kl + 4} U {kl + a^} = SpecVn U SpecV22, 

/ci , /c2 G R, whereas the spectral density over a large interval of length L is : 

for fc| > 4 - (T^ 
for A;| < 4 - cr^ 

From these data, the heat and generalized zeta functions can be readily computed : 
Tre-^^ = Tre-^^" + Tre"^^- = ( e"^^ + e"'^'^) quad, 

• Generalized zeta function of K, 

In the case of the one-component topological kink, computation of the heat and generalized 
zeta functions is easy because K is diagonal. Moreover, Ku and ^22 are respectively the Hessian 
operators for the A(0)^ kink and the sine-Gordon soliton (in the second case shifted by — 1). 
Therefore, we shall take advantage from the work performed in [Q. The spectrum of the Hessian 
operator for the TKl kink is : 

Spec£ = {0, 3, a' - 1} U {kl + A}k,e^ U {kl + a^}k,^^ , 

and the spectral density and the phase shifts for the continuous spectrum read : 

p^ik,, k,) = p,(fc„ k,) + - + J 

3/^1 1 
^i(^i) = — 2arctan- , 52(^2) = 2arctan— 

respectively. Besides the continuous spectrum there are three bound states with eigenvalues 0, 3 
and (7^ — 1. The eigenfunction of eigenvalue is the translational mode. There is a second bound 
state in the ei direction, but the third bound state points along the 62 axis. Note that to develop 
a non-zero 02 component is energetically favorable if cr^ < 1 and this process costs energy when 
> 1 and the TKl kink is stable. The second zero mode that occurs at cr^ = 1 is the signal of 
this phase transition. 



Pv{h,k2 



mL 

mL 
2\/2- 



TT 
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The heat function 



Tre"'^^'^ — 'jj.g-/3-P/Cii _|_ rj^g-/3K:22 
dki 



mL_ ^ J_d6ik^\ ^_^;.2 
y/SiT 27r rf/ci y 

mL 1 d6{k2)\ 



+ 



v/Svr 27r (iA;2 / 



is therefore equal to : 

Tre-^^^ = Tie-^^ + fe-^^ + e'^'^'-i)^) Erf [v^] - Erfcpv^] 



where Erf and Erfc, are respectively the error and complementary error functions. The Mellin 
transform provides the corresponding zeta function: 



Cpds) = Cv{s) + 



3.+ 



2 IP n 1 3 n 1 2 



r2-Fi[|, -5 + 2' 1' P"] 



r(.) 



in terms of Gauss hypergeometric functions of the form 2-^1 [o, c, d], |jT6|. 
Applying these results we obtain : 



K 



— lim 



2/i' 



2\ « 



X 



-2Fl[|,S + 



- ^ ^1 



r(s) 



which is still a divergent quantity. The mass renormalization counter-terms add another divergent 
quantity : 

,2\ 



A K V 2n /2/i- 

Ase = - lim hm — [ 



r(s + i) 
r(.) 



[3Cv,,{s + I) + Cvjs + 1)] 



(9) 



and we obtain: 

Aie^ + A2e^ = -hlim 



+ — hm 

2 £^0 



my2 /2/i 



£^0 

2N- / 1 




2/^2 



271 \ 



i-t IL2' ^' 2' 3J 



3 1 
— + 



2\e 



_^2^ie2-^l[2'^' 2' <72- 



l ^ 3 1 



It is appropriate to consider the contributions from /Cn and /C22 separately before taking this 
limit: 



From /Cii we obtain 



+ 



3hm 
- hm — 1= h 

e^o 2V27re 

hm 



2v^7r 



2/i2 



-4 + 37£ + log(3 ■ 2^) - 3 log ^ + 3V^(|) - 2^^ [i 



2' 2' 3J 
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and 



lim — 1= — 

2v/27re 



hm 



2V2' 



-37£; + 31og^-61og2- 37/^(1) ) + o(e) 



Here, 7^ is the Euler gamma constant; ip is the Psi(Digamma) function, and is the 
derivative of the Gauss hypergeometric function with respect to the second argument, |TB 



From /C22 we obtain 



and 



A 



+ 



hm 

— hm — — — 

hm 



2^2: 



TT 



+ 



, 2/i2 

log + ^og{a 
m'^ 



1 



' 2' 1-0-2 



0(e) 



hm — ^ — 

^^0 2v/27re 



+ 



hm 



Gathering all this together, we finally find: 



--iE + log 



2/^2 



m^ 



loga^-V^d) +o{e 



Ai£^ + A^e 



K 



A^pe^ + A, 
hm 



e 



+ Af + A 



hm 



4 + ln^ + 2F{[i,0,f,-i]+ln(- ^ 



vr 



(3 



TT 



12 



+ (1- 



1 arcsin— ) 
a 



(10) 



Note that AMtki(o'^) = Aie^ + A^e^ acquires an imaginary part if cr^ < 1 and in this latter 
regime the quantum TKl kink becomes a resonance. Formula (|T0|) shows one interesting pattern 
in the one-loop correction to the kink masses: First, the arguments of the logarithms are always 
quotients of the eigenvalue where the continuous spectrum starts by the energy of the bound 
eigen-state. Second, the fourth argument of the derivatives of the Gauss hypergeometric functions 
is always minus the inverse of the bound -state eigenvalues. 



3.2 The quantum topological kinks: asymptotic expansions 

The difficulty in computing the quantum correction to the mass of the TK2 kink lies in the fact 
that the 7i Hessian operator at this configuration is a Schrodinger operator for which the 2x2 
matrix- valued potential V[x) is non-diagonal. Therefore, the spectral problem of li, is very difficult 
to solve and we are led to use asymptotic methods to calculate the quantum TK2 kink mass. 
In the MSTB model there are two regimes: 

1. 0"^ < 1: The TK2 kink exists and is stable. In this case we have that 

- , , ^-r , . ^_tanhcrx 

is an odd function of x such that [ai]i2('^^) = [ai]2i('^) = 0. Bearing in mind that V\\{x) = 

~ clsh4i ^"^^ ^22(3;) = - "^coswZJ ' present the coefficients [a„]aa(H) up to no = 11 with 
the help of Mathematica in the Tables 1 and 2 for several values of a: 
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(7 = 


0.5 


a = 


0.7 


(7 = 


0.8 


n 


KliifH) 


K]22('H) 


KliiCH) 


[an]22(W) 


KliifH) 


K]22(W) 


1 


17.9984 


-9.99909 


14.2285 


-3.02857 


13.2000 


-0.4000 


2 


34.9975 


16.3313 


27.5051 


4.95579 


25.6320 


2.42133 


3 


47.9636 


-16.8309 


38.5365 


-3.88117 


36.3858 


-1.46193 


4 


49.4219 


13.2491 


40.9039 


2.58248 


39.2877 


0.984184 


5 


40.6415 


-8.41779 


34.5847 


-1.43041 


33.7529 


-0.53935 


6 


27.7911 


4.4983 


24.2279 


0.696572 


23.9672 


0.267282 


7 


16.2632 


-2.08006 


14.4768 


-0.30375 


14.4819 


-0.117869 


8 


8.31759 


-0.849979 


7.53962 


0.119982 


7.61204 


0.0464771 


9 


3.77771 


-0.295903 


3.47957 


-0.0425253 


3.53989 


-0.0162638 


10 


1.51992 


0.0884623 


1.44002 


0.0132872 


1.47558 


0.00497705 



Table 1: Heat function coefficients for several values of the deformation parameter such that a < 0.9 . 





(7 = 


0.9 


a = 


0.95 






0- > 1 


n 


K]ii(w) 


K]22CH) 


KiiiCH) 


K]22(W) 




n 


[a„]ii{IC) 


[a„]22ilC) 


1 


12.4889 


1.91111 


12.2211 


2.97895 




1 


12.0000 


4.0000 


2 


24.5218 


1.67378 


24.1951 


1.952444 




2 


24.0000 


2.66667 


3 


35.3368 


-0.241656 


35.16.16 


0.310882 




3 


35.2000 


1.066667 


4 


38.8216 


0.341514 


38.9553 


0.230323 




4 


39.3143 


0.304762 


5 


33.8715 


-0.16798 


34.2227 


-0.053986 




5 


34.7429 


0.0677249 


6 


24.3610 


0.0863412 


24.7475 


0.0372813 




6 


25.2306 


0.0123136 


7 


14.8742 


-0.0370868 


15.1758 


-0.0145382 




7 


15.5208 


0.0018944 


8 


7.88516 


0.0140749 


8.07351 


0.00547229 




8 


8.27702 


0.000252587 


9 


3.69263 


-0.0046464 


3.79191 


-0.00173035 




9 


3.89498 


2.9713410-5 


10 


1.54847 


0.001313 


1.594006 


0.0003535 




10 


1.63998 


3.12591 10"^ 



Table 2: Heat function coefficients for several values such that a > 0.9. If a > 1 the coefficients coming 
from K. are shown . 
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We also need to know that J = 1: there is only one translational zero mode in the spectrum 
of Ti, and that vl = A , V2 = u"^ ■ Plugging this information in the asymptotic formula 
we obtain the values shown in Table 3 for AMtk2- 

2. cr > 1: The only solitary wave is the TKl kink, which is stable. To determine the quantum 
TKl kink mass by using the asymptotic method, the general formulas must be applied to 
IC. Now, Vi2{x) = V2i{x) = and trivially [ai]i2(^) = [an]2i(^) = 0. The [a„]aa(^) must be 
read from the formulas in the Appendix applied to Vu{x) = ^ ^^^^i^^ and ^22(2;) = — ^^^^a^^ , 
and are shown in Table 2 up to = 11. 



a 


AMtK2 


0.4 


-1.103270;im 


0.5 


-0.852622;im 


0.6 


-0.68900i;im 


0.7 


-0.583835;im 


0.8 


-0.524363;im 


0.9 


-0.505708/1771 


0.95 


-0.511638/7777 



a 


AMtki 


1.0 


-0.528311^777 


1.2 


-0.518426/7777 


1.4 


-0.509645/i777 


1.6 


-0.502291^777 


1.8 


-0.496369^777 


2.0 


-0.49172^777 


2.5 


-0.484183/7777 


3.0 


-0.480101/i777 


4.0 


-0.476181/i777 



Table 3: One-loop corrections to the mass of topological kinks 

Note that in this case the coefficients of the asymptotic expansion of the generalized zeta 
function of K, are independent of a. Nevertheless, the quantum correction to the TKl kink 

mass depends on a through the factors ^^"a^-'a^^ . Again j = 1: there is only one translational 

'"2 

zero mode in the spectrum of /C, and vj = A , V2 = cr'^ . Application of formula (|^) provides 
the values of AMtki, also shown in Table 3 for several values of a. 

We summarize the results obtained in this Section in the next Figure, where AMk is depicted 
as a function of a. It is understood that for a > 1 we plot AMtki and for < a < 1 AMtk2 
is represented. The continuous line shows the exact value of AMtki as a function of a, whereas 
the dots correspond to the answer provided by the asymptotic method to AMtki, (whites ), and 
AMtk2, (blacks ), for several values of a. 
We observe the following facts: 

• In the a > 1.2 range, the approximation established by the asymptotic method is extremely 
good : the discrepancy with the exact answer is compatible with zero. 

• In the 1 < 0" < 1.2 interval the error is larger and can be estimated exactly in terms of 
-S/c*(~i), ^no,/c(~|), Bv{~j) and i?v^(|). The closer the value of cr to 1, the smaller the 
first eigenvalue in the spectrum of }C* and the larger i?£*(— |). The asymptotic method is 
better when fiuctuations in the 02 direction of the TKl kink cost more energy. 
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-0 . 4 
-0 . 5 
-0 . 6 
-0 . 7 
-0 . 8 



AMk 



-e-Q- 



q 

4 



Figure 1: Quantum correction to the mass in the MSTB model: TK2 (t) < cr^ < 1, hlack dots) and TKl 
(a^ > 1, white dots). 



The point o" = 1, where the error is maximum, is singular. There is a second zero eigenvalue 
which is the signal for an instability-type phase transition on the TKl kink. The co-existing 
unstable TKl and stable TK2 kinks for a < 1 coalesce to a single TKl kink at o" = 1, which 
becomes stable for a > 1. 

In the < a < 1 regime there is no way of estimating the error in the approximation to 
AMtk2 obtained by means of asymptotic methods because the exact value is not known. 
We expect, however, that the result will be better for smaller values of a because of the 
same qualitative argument as above: the cost in energy for climbing from the TK2 to the 
TKl kink is larger with smaller a. 

The value cx^ = 4 is a very special one. In this case the system is A/" = 2 pre-supersymmetric 
in the sense that there exists a super-potential : 



>i,(/>2je2 



where 



2 

2 ) 



n,(P2 



^102 + — 



and 



If + 202 



1 dW^ dW^ 1 dW^ dW^ 



2 (901 (901 ^ 2 (902 902 2 (901 90i ' 2 902 902 
This is no more than the super-potential of the holomorphic Wess-Zumino model. 



1 dW^ dW^ 1 dW^ dW^ 
\ 

2 



see 



r7[], and we have thus calculated the quantum correction to the kink in this important 



system. 



4 Quantum kinks in a deformed 0(N) linear sigma model 



Let us finally consider an 0{N), N > 3 model determined by the potential energy density: 



m 



N n 

a'' 



(11) 



1=2 
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which generahzes the system discussed in Sect. 3. Without loss of generality, the deformation 
parameters are chosen in such a way that 02 < as < ... < 01 n- Passing to non-dimensional 
variables, the potential is: 



1 N 2 



1=2 

\N . 



with a2 < 0-3 < . . . < gm- The "internal" symmetry group is (Z2) , the vacuum classical 
configurations are : 



(t)v{x,t) = ±ei , ipv{y,y^) 



m 



ei 



and the vacuum orbit and the vacuum moduli space are respectively: Ai = -^fy = Z2, Ai = ^ = 
point. 

There is a very rich manifold of kinks in this model, see We shall only study the one- 
component and two-component topological kinks, which are loop kinks and candidate stable soli- 
tary waves in this system: 

• Topological kinks and antikinks with one non-null component: 

0TKi(a;,t) = ±(tanhx)ei , ?/'TKi(y, 2/°) = ±-^(tanh--^)ei . 

• Topological kinks and anti-kinks with two non-null components. For each n from 2 to 
such that cr„ < 1, we have : 

</'TK2f7„(a;,i) = ±[(tanhcr„2;)ei ± a„(sechcr„x)e„], 
V^TK2a„(l/,/) = ^[m(tanh^)ei ± v^(m2^"a|)(sech^)e;] 



where cr„ = a/1 — cr^. 

The kink and vacuum solutions have classical energies: E['?/vpKi] = E['?/>rK2cr,J = 2(t„(1 — 

23—* — * — * — * 

^)^, E[?/'y] = 0. Thus, E[^TKi] > E[V'TK2aiv] > • • • > E[V^TK2a2]- The lowcr bound in energy in 
the topological sector of the configuration space is attained by the TK2cr2 kink. If (12 > 1, only 
the TKl kink exists and is stable. If cTj < 1 < cTj+i, the kinks TKl, TK2cr2, . . ., TK2aj exist, but 



only the TK2(T2 is stable, see [IS 



The Hessian operators in non-dimensional variables are as follows: 
• Vacuum: 

V = dmg{Vu) 



Vii = -— + 4, Vu = -— + af, I > 2 



TKl: 



/C = diag(/C;,) 

6 ^ rf2 ^ 2 , 

/Cii = -3-^ + 4 K,ii = -— + ai l>2 

dx-^ cosh X dx^ cosh X 
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TK2(T„: 



(__£_ _|_ ^ _ 2(2+o-„) tanhgna 
cosh^ anX " cosher,! x 

" cosho-„a; dx^ ~ n~ cosh^ crno; 

where the hat over an index means that this index is absent and rows and columns have 
been rearranged in such a way that 7ij^]j acts on the small deformations of 0i, 0„ and 7i^^ 
on the remaining ones. 

These operators are direct sums of the Hessian operators used in Sect. 3. We can therefore take 
advantage of the calculations already made to give the quantum corrections to the kinks of the 
present model. 



4.1 The quantum TKl kink: exact computation of the semi- classical 
mass 

The relevant quantities are: 

• Generalized zeta function of V: 



N 



N 



Cv(5)=Cvn(^)+ECva(^) 



Vis 



1=2 



mL I 1 1 



Generalized zeta function of PIC : 



1 ^ 



2 1 1 ' 

r2i^i[|,s + -,|,-|] - — 



3^+2 

9 



V{s) 



r{s) 



Applying these results, we obtain : 



2 1 ^\m^ 



—2Fi[l,s + -,l,-l] - ^ 



h 1 f2u' 

+ - hm — -h- 
2 ^/^^ \m? 



N 



1=2 



3"+2 



2 1 

r2-^i[|,s + -,|,4 



r{s) 



-afY ' 2 



The mass renormalization counter-terms are : 



N 



. K , , 2n /2u2V+2 r(s + 1)^ ^ ^ 

A2£^ = - lim. lim -[A,] 4^[3Cvn(5 + 1) + E Cv.(^ + !)]• 



r(.) 



/=2 
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Therefore, 



2V27r 
hm 



A / 2 \ 

4 + In - + 2F[[l 0, + (^^(^) + 0' i' ^] J 



;=2 

N 



7rV2 ^ 

The TKl kink is a bona fide quantum state only if erf > 1, VZ > 2. 



(12) 



4.2 Quantum masses of two-component topological kinks: asymptotic 
expansion 

In the cr| < 1 regime only the TK2(72 kink is stable. In this sub-section, therefore, we shall present 
the calculation of the one-loop quantum correction to the TK2(72 kink mass. Since Tre""'^^^^^ = 
'jYg-/3Wjv£, _|_ r^^-PHjy ^ encounter two old friends: the heat functions arising respectively in 
connection with the TK2 kink in the MSTB model and the soliton of the sine-Gordon model. 
Using the information collected in previous Sections, we find: 

• Casimir energy : 



hm 

-|- — ;= lim 
2V2 



s- 



Mass renormalization energy : 



=3 

We finally obtain the answer: 



* \ * a=l n=l 



^2n-2 



or 



*^ \ ^ a=l n=l " 
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In both formulas the contribution coming from 7i)vi) '^^^ t)e read from the information on the 
quantum correction to the TK2 kink mass in the MSTB model, collected in the Tables and the 
Figure of sub-Section §.3.2 . A similar formula would show us that AMtk2ctp V/ > 2, receives an 
imaginary contribution any TK2ai quantum kink state is therefore a resonance. 



5 Further comments 

The family of deformations that we have treated admits a A/" = 1 super-symmetric extension. The 
super-potential W{(t>) in the N = 2 case is: 

W{ct>) = ±7 (0? ± ay + 0i [^(0? + ^It a0i + a') - 1] . 

For the special value a = 2, this system also admits the super-potential mentioned at the end 
of sub-Section §3.2 because it becomes the Wess-Zumino model. If iV > 3, it is very difficult to 
write the super-potential in Cartesian coordinates in the internal space; nevertheless, passing 
to elliptic coordinates one obtains easy expressions for the super-potential, see [FHl - 
The super-symmetric extensions also include Majorana spinor fields: 

x{x^) = ^ X2{x'^) ) ' ~ Xa ; a = 1,2 

Choosing the Majorana representation 7° = cr^, 7^ = ia^, 7^ = of the Clifford algebra {7^^, 7''} = 
2g'^'^ and defining the Majorana adjoint x = X*7°) the action of the super-symmetric model is: 

S = dx^ [d^$d^'$+ iXl^d^X - VWVW - flwx] , 



a=\ a=l b=l 



The M = 1 super-symmetry transformation is generated on the space of classical configurations 
by the Hamiltonian spinor function : 

Q = y"da;|7Vx5M'?+^7°xW} . 

The components of the Majorana spinorial charge Q close the super-symmetry algebra : 

{Qa, Q/3} = 2{Y^XpP, - '^^llpT. (13) 



Their (anti)-Poisson bracket is given in (|T3D in terms of the momentum and the topological 
central charge T = \ j dW\ = \ j VWd(^\ . 

The chiral projections Q± = ^^^Q and x± = ^^^X provide a very special combination of the 
super-symmetric charges: 



Q^ + Q_= I dx{ (x+ - X~)^ - (X- + X+)^W 
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+ Q- is zero for the classical configurations that satisfy ^ = =FVVr and x± = which are 
thus classical BPS states. In Appendix B it is proved that the stable TK2 kinks are such BPS 
states and besides the small bosonic fluctuations one must take into account the small fermionic 
fluctuations around the kink in order to compute the quantum correction to the kink mass in the 
extended super-symmetric system. The fermionic fluctuations around the kink configuration lead 
to other solutions of the field equations if the N x N matrix Dirac equation 

is satisfied. We multiply this equation for the adjoint of the Dirac operator : 
{-zYd, + ^W{$k)} + Kw{^k)} ^FXix.t) = , 

and, due to the time- independence of the kink background, look for solutions of the form: Spxi^) t) = 
u;)e*^*. This is tantamount to solving the spectral problem 

- + lwi$K)^Wi$K) T z7'Viy (0k) V ® Kwi$K) I frix; cu) = u^pix- u) . 

Projecting onto the eigen-spinors of ^7^, 

we end with the spectral problem: 




for the same Schrodinger operator that governs the bosonic fluctuations. 

Therefore, generalized zeta function methods can also be used in super-symmetric models to 
compute the quantum corrections to the mass of BPS kinks. Great care, however, is needed in 
choosing the boundary conditions on the fermionic fluctuations without spoiling super-symmetry. 
We look forward to extending this research in this direction. 



1 / fp^ixiuj) - fpjx]uj) \ 

2 I -fHx^^^) + fpix^uj) I 



A Appendix: the matrix heat kernel expansion 

In this Appendix we describe how to find the coefficients of the asymptotic expansion of the heat 
kernel associated with K, and show the explicit expressions for them up to third order; for a more 
complete treatment see |T^. Formula (|^) in the text tells us that the matrix elements of A{x, x'; (3) 



satisfy the coupled PDE: 

N 



d_ ^ ~ ^} [^Uix,x';/3) = J2yacix)[AU{x,x';/3) + (v^ - vl)[AUix, x'; p) 
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with the initial condition: Aab{x,x';0) = 6ab- Plugging the power expansion of Aab{x,x'; P) into 
this system of equations we find the recurrence relations: 



n[an]ab{x,x') + {X - X')^^^^^{X,X') = ^ ) +y^^Vac{x))[an^l]cb{x,x') 

c=l 

+ {v"^ - vl)[an-l]ab{x,x') . 

In order to take the x' — ^ x limit properly, we introduce the notation: 

")l^„],.(.)=lim%k(,,,') . 

Then, the recurrence relations in the x' = x limit become: 

if ^ - 

[an]ab{x, X) = - < [An-l]ab{x) + ^ Vac{x) [an^l]cb{x, x) + {v^ - vl) [an-l]ab{x, x) 



N 



n 



c=l 



We also need the secondary recurrence relations among the ^^''[A^ab{x) derived from /c-times 
differentiation of the primary recurrence relations above: 



^'\AX,{X) = ^ Y'^'\A^-Aab{x) ^j^i^i]) '\An-lUx) + {vl - vlt\A^-Aab{x) 

I c=l j=0 ' 



Notice that '^'^^A^Uix) = lim,,^, = S'^^ab- Thus, the ^""^AnUix) and, hence, the [an]ab{x,x) 

can be generated recursively . The three first coefficients are : 

[«2]„5 (^) = + \ [V\b (^) + - vl)Vab{x) 

N„. (^) = + ^ [V^'\x)V{x)]^^ + ^ [V{x)V^'\x)]^^ + ^ [y«(x)y«(a:)]^, 

1 r.-.,i . . 1 



+ I [y\b (-) + Y2^-i - 1) {vS:\^) + 2 [v\, {X)} 

1 1 ^ 



c=l 



We mention that, as happens in the scalar case [^, the diagonal terms [an]aa{x,x) can be 
interpreted as the densities giving the infinite conserved charges of a matrix Korteweg-de Vries 
equation; namely : 

dV - dV dV - dW 

dt ~^ '^^^ dx ~^ dx^'^ ~^ dx'^ ^' ^ ^ 

where now the matrix potential evolves in "time" t, V = V{x,t). The reason is that can be 
written Lax equation 

Li + [L,M] =0 
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for the operators 



M = A—-3V—-3—V + Bit) . (16) 
ox-* ox ox 



with i?(t) arbitrary. Therefore, standard arguments guarantee that the time evolution ruled 



by ([MD produces an uniparametric isospectral transformation of the Schodinger operator (p!5|). 
Because the integrals [a„]aa are determined by the spectrum of (|15|), their invariance follows. 

B Appendix: BPS and non-BPS kinks 

This Appendix is devoted to characterize respectively the TKl and TK2 kinks as non-BPS and 



BPS states in the N = 2 case, see Reference [181. We have seen in Section 55 that the model 



admits four super-potentials. If /3a = 0, 1 , a = 1, 2 , the four super-potentials 

W^^^'^-\^,, 02) = (-l)^V(0i + (-l)^^^)' + 02[^(0? + 4>l- (-l)^V0i + a') - 1] , (17) 



satisfy: 

1 A aiy^^i''^^) 1,,2 .2 .^2 

-> 7. = -(0T + 09-l) H 

2^ d(t)a d(P„ 2^^^ ^2 ^ 2 

The energy for static configurations reads: 



1 /■ ^(d(j)a dW'^'^''(^^^\ fd(f)a r J^dclyadW^'^'''^^^ 

E-- dx2_^{— — — — + / dxl^ 



- o \ dx d(j)a J \dx d(j)a J J dx 

The BPS kinks are the solutions of the ODE first-order system : 

_ + 02) (^^ ^ (-l)fea )+0,02] 

^2'+f01 + (-l)''V^' 



(i0i 


5P^(/3i,/32) 


(ix 


901 


d(j)2 


5|y(/3i,/32) 


dx 


(902 



5H^(/3i,/32) 02 [_^2 ^ 0^2 ^ ^^^^^ ^ {-If'cr)] 



^02'+(01 + (-l)^V)' 

From the values of the super-potential at the vacuum points 

^(/3i,/32)^±l^ 0) = {-if '{I ± (-1)^V)(-1 + ^(1 T {-If'a + a^)) 
we calculate the Bogomolny bound 

E^P' = I j rfW^(A,/32)| ^ |i^(/3i,/32)(i^o) - 1^(^^'^^)(-1,0)| = 2a (^1 - 



(19) 



(20) 
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which is saturated by the solutions of (0,|T^), the BPS kinks. 

In the derivation of (pO]) we have used Stokes's theorem. The foci (0f = —{—l)^^a, 0f = 0) of 
the elhpse (pf + = 1 in the "internal" space are branching points of iy(^i'^2). Therefore, 
we are legitimated to use Stokes's theorem - and the Bogomolny bound is reached- only if the 
kink configuration does not cross any of the two foci above mentioned. The TK2 kinks live on 
the semi-ellipses = ±a^^/l — (0™^)^ and everything is fine: they are BPS kinks. The TKl 
solutions are more involved. If 02 = (|^) reduces to: 

-ifni-0D |l^!i~!|';!i = -(-If ni-0?) (e(0i + ai-lf^) - e(-0, - a(-lf ^)) 



where Q{u) is the Heaviside step function. Thus, 0™^(x) = tanhx is not solution of the first-order 
equations on the whole real line; if 0i + (— l)^^cr < 0, we must choose Pi = 0, and Pi = I otherwise. 
On the half-line x E {—oo, — arctanh[(— l)^2cr]] the TKl solution is the flow line of gradVT^^'^^-' but 
it becomes the flow line of giadW^^''^^^ on x G [— arctanh[(— 1)^2(t], oo). One can easily check that 
the TKl kink is a proper solution of the second-order equations with energy given by a piece-wise 
application of Stokes's theorem: 

3 

In a super-symmetric extension of this model the corresponding state would be not annihilated by 
any combination of the super-symmetry generators built from one of the ^/(^I'^a) super-potentials 
and, therefore, the TKl kink is a non-BPS state in this system if cr^ < 1. 

It is shown in Reference |jl8| that only the stable TK2 kinks are BPS kinks in the > 3 



models; the proof has been performed using elliptic coordinates for the field variables. 
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